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Ah St nact 

The ai r< rriO parameter estimation problem is 
iS ed to illustrate the utility of parameter esti- 
mation, which applies to many engineering and 
sc ienti tic. fields. Maximum likelihood estimation 
r« ( »s been used to extrac: stability and control 
dprivatives from flight data for many years. This 
oauer presents some of the basic concepts of air- 
craft parameter estimation and briefly surveys 
the literature in the field. The maximum likeli- 
hood estimator is discussed, and the basic con- 
cepts of minimization and estimation are examined 
for a simple simulated aircraft example. The cost 
functions that are to be minimized during estima- 
ti 0 n are defined and discussed. Graphic represen- 
tations of the cost functions are given to illus- 
trate the minimization process. Finally, the 
basic concepts are generalized, and estimation 
from flight data is discussed. Some of the major 
conclusions for the simulated example are also 
developed for the analysis of flight data from the 
F-14, Highly maneuverable aircraft technology 
( HiMAT ) , and space shuttle vehicles. 

Nonenc 1 ature 

A, H, r , , system matrices 

n, F , G 

a y lateral acceleration, g 

; ^ coefficient of rolling moment 

r fn efficient of pitching moment 

; coefficient of yawing moment 


L rolling moment, divided by ! x , 

deg/sec 2 ; or* iteration 
number 

L ' rolling moment, ft-lb 

} yj rolling moment due to yaw yet, 

ft-lb per jet 

M Mach number 

N number of time points or cases 

n state noise vector; or, number of 

unknowns 

p roll rate, deg/sec 

q dynamic pressure, lb/ft 2 

Re Reynolds number 

r yaw rate, deg/sec 

t time, sec 

u control input vector 

V total velocity, ft/sec 

x state vector 

z observation vector 

z predicted Kalman-filtered estimate 

a angle of attack, dey 


^dyn 

fy 

t { - ) * y(-) 

G(i* 


< ■ *y * 

1 z ♦ J x7 


equivalent dynamic directional 
stabi 1 i ty 

coefficient of sideforce 

general functions 

measurement noise covariance matrix 

approximation to the information 
matrix 

moment of inertia about subscripted 
axis, slug-ft 2 

cost function 


*Substant ial portions of this paper are 
taken from two publications of the author. 
Refs. 1 and ?. 

**Chief, Fluid and Flight Mechanics Branch. 
A I AA Fellow. 
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<$DE 

«r 

n 

M 

a 

T 


angle of sidesl ip, deg 
time sample interval, sec 
control deflection, deg 
aileron deflection, deg 
differential elevon deflection, deg 
rudder deflection, deg 
measurement noise vector 
mean 

vector of unknowns 
standard deviation 
time, sec 


Mm paper i* declared a work of the U S 
i.ovrrnmem and therefore it >n the public domain. 


♦ 

if 

Subscript s 


transition matrix; or, bank angle, 
deg 

integral of transition matrix 


are time functions. Autoregressive moving average 
(ARMA) models are also widely known; however, 
discrete-time ARMA models can readily be rewritten 
as linear state-space models,* so the discussion 
of state-space models presented in this paper is 
applicable to ARMA models. 


P, r , a, 0, 

$nE. V 


partial derivative with respect to 
subscripted variable 

bias ; or, at time zero 

minimum value 


Supersc ript s 

predicted estimate 
estimate 

* transpose 

Introduction 


It is difficult to present a topic as special- 
ized as aircraft parameter estimation in a way 
that will interest a generalized audience of 
mathematicians , scientists, and engineers. The 
approach here is to portray parameter estimation 
as a specialized "curve-fitting” technique that 
can be applied to a broad class of problems. Much 
effort is expended in a variety of disciplines on 
a form of curve fitting, more specifically, the 
correlation of observed or inferred data with an 
assumed (though perhaps in a high- or infinite- 
dimensional space) mathematical model that is 
based on phenomenological considerations. This 
broad class of problems is referred to as system 
identi f ication. 

[he application of system identification, 
sometimes referred to as the inverse problem 
(paraphrased as. Given the answer, what is the 
question?), presumably goes back to prehistoric 
times as humanity tried to master the environment 
by understanding, based on observation, certain 
phenomena (probably simple ones). Many of the 
physical laws stated by the Chinese, Egyptians, 
and Greeks were based on the same principles as 
are currently used in system identification. 

Through advancing technology and mathematical 
rigor, we can apply much more sophisticated tech- 
niques for making observations and for deducing 
the underlying phenomenology, but the basic prob- 
lem of system identification remains the same. 

For most physical systems, information about 
the general form of the system to be identified 
often can be derived from knowledge of the system. 
The most widely applied subfield of system identi- 
fication is parameter identification, where the 
form of the mathematical model is assumed to be 
known. The model (an explicit function, a polyno- 
mial expansion, a look-up table, a finite-state 
machine defined for application of artificial 
intelligence, or many other forms) contains a 
finite number of parameters, the values of which 
need to be deduced or identified from the obser- 
vations. One of the favored forms of the model 
for the most successful application is the state- 
space form (a rigorous treatment of state-space 
forms is given in Ref. 3). State-space models are 
very useful for dynamic systems, in which responses 


An assumed model wi 1 1 not be an exact repre- 
sentation of the system no matter how carefully 
its f ' nn is selected. The experimental data will 
not be consistent with the assumed model for any 
parameter values. The model may be close but will 
not be exact, if only because the measurements or 
observations will be made with real, thus imper- 
fect, sensors. Errors in observations and in the 
model need to be evaluated in determining the 
unknown parameters of the model. So tne objective 
becomes the application of the "best' 1 model (in 
some sense), instead of the correct model, to find 
the "best" estimates Oi the unknown parameters; 
this process is referred to as parameter estima- 
tion. The currently favored approach to parameter 
estimation, and the one discussed in this paper, 
is to minimize the error, in the least squares 
sense, between the model response and the actual 
measured response; the estimates resulting in the 
minimum error are the "best" estimates. The 
theoretical formulation^ and application^ of the 
output error technique (which is a maximum likeli- 
hood technique that is used throughout this paper) 
have been thoroughly documented. 

Although the applications described in this 
paper pertain to aircraft, the techniques have 
been successfully applied in other fields where 
the mathematical model and observations are ade- 
quate. Parameter estimation may sound like one 
more arcane subject, but it has application in any 
field where observations must be made to agree 
with the assumed physics of a problem. There are 
many obvious applications in a variety of fields, 
such as, spacecraft dynamics, gravitational per- 
turbations, fluid dynamics and mechanics, optimal 
control, and guidance. 

The application of the maximum likelihood 
technique for parameter estimation of aircraft 
coefficients demonstrates a successful application 
of system identification technology. Analysts 
in the aircraft community accept and use system 
identification techniques on a routine basis. 
Although there are isolated problems (primarily in 
extending the application to more difficult flight 
regimes, such as where the aircraft is dominated 
by poorly understood separated flow), there f s 
little doubt that the basic application is highly 
successful. Contributing to this success are a 
wel 1 -understood, time-tested, 7 *® physically 
derived model form that is reasonably represent- 
ative of the true vehicle in most flight regimes; 
high-quality measurements of several relevant 
states; the ability to apply Inputs specifically 
for system identification; and engineers familiar 
with system identification, aerodynamics, aircraft 
equations of motion, and the associated aerodynamic 
coefficients. 

This paper first presents a brief survey of 
the contributions to system identification, and 
specifically aircraft parameter estimation, up to 
1980, when the maximum likelihood technique began 
to completely dominate the field. (Refs. 6 and 9 
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give a broad view of contributions since 1980. 

Ref. 9 is a bibliography of nearly 500 books, 
papers, and reports related to parameter esti- 
mation.) Some common uses of the estimated 
parameters are then discussed. The technique 
used tor parameter estimation is then described, 
followed by an examination of the computat ional 
details and cost tunctions involved in error mini- 
mization. Mnally, applications of the technique 
for improving hi gh-performance aircraft and the 
space shuttle are described. 

History o f P a r am eter Es t iroat i o n 1 0 1 980 


Pe n era I Sy s t em_ Men ti fication 

Thf* transition from hit-and-miss, rule-of-thumb 
system identification to mat hemat i ca I I y sound 
approaches has been gradual; certainly no single 
germinal work can be referenced. Gauss, 10 in 
1H()9, discussed the inverse (system identifi- 
cation) problem and suggested some statistical 
approaches that are relevant even today. The 
discussions of Douglas, 11 in 1940, and Gelfand 
and Lpvitan,^ in 1951, pertaining to the inverse 
problem certainly qualify as truly significant 
works contributing to the state of the art. The 
formulation by Feldbaum, 13 one of the more sig- 
nificant works aimed at the current direction of 
invest igation, is somewhat different than others 
discussed, but he did look at identi fication and 
control of the system as a single problem, the 
“dual control" problem. During the 1960s, a 
plethora of publications was evidence of increased 
interest in problems of this type. Much of this 
interest was stimulated by the well-known early 
works of Kalman. 

The bulk of general system identification 
theory and application up to 1980 has been sum- 
marized in several excellent survey papers. 14-1 7 

The system ident i f i cation problem can be 
divided into two major subsets: deterministic 

(without state noise) and nondeterministi c (with 
state noise). There are two classes of tech- 
niques for identi fication of nondetermi m st ic 
systems: the Kalman filter (or more generally, 

the extended Kalman filter) technique and the 
maximum likelihood technique. Many precise 
applications do not truly fall into these classes, 
but they do tend to mimic one of the two tech- 
niques. The extended Kalman filter (discussed 

o 

by Ast.rom 1 ^ and Kashyap^ 9 ) has been widely 
applied; however, this paper primarily examines 
the maximum likelihood estimator, proposed by 

8a lakrishnan 20-22 an d developed in Refs. 23 and 24. 
Ai rcra ft I dentification 

In the following chronological survey of 
investigations that led to the development and 
widespread acceptance of the maximum likelihood 
estimation technique for aircraft coefficient 
estimation, the more st raight forward deterministic 
analysis is discussed first, followed by a brief 
discussion of nondetermi ni st i c analysis. Some of 
the investigations in estimation of unknown coef- 
ficients from aircraft dynamic response data are 
contained in Refs. 25 and 2b. The National 


Advisory Committee on Aeronautics (NACA) had been 
publishing reports on stability derivatives 
(coefficients of the differential equations of 
motion) since the early 1920s. (The reports by 
Norton?* 8 involved the identification of frequency 
and damping ratios from flight data.) 

Deterministic Analysis . The sophistication 
and complexity of the methods used to estimate 
unknown coefficients from aircraft dynamic flight 
responses have increased over the past 40 years. 

In the late 1940s and early 1950s the frequency 
response methods (including steady-state oscilla- 
tor analysis 27 and Fourier analysis 2 ^) increased 
in popularity in aircr^t analysis and in other 
applications, ihese methods yield the frequency 
response of the vehicle but not the coefficients 
of the differential equations. Attempts were made 
to extract tnese coefficients by selecting values 
of the aircraft coefficients that resulted in the 
best fit of the frequency response. 29 * 30 Regres- 
sion techniques, such as linear least squares 31 
and weighted least squares 30 techniques, were 
also applied to flight data at about that time. 
Unfortunately, regression techniques give poor 
results in the presence of measurement noise and 
yield biased estimates. The time vector tech- 
nique 32 has also been applied to flight data; 
however, it yields an incomplete set of coef- 
ficients, and the types of responses that can be 
analyzed are restricted to fairly simple motions. 
Analog matching techniques 32 * 33 (time consuming 
and somewhat tedious) have also been applied to 
flight data but are limited because resulting 
estimates vary with the skill and technique of 
the operator. Comparisons of these early tech- 
niques 34 * 36 showed that a more complete method 
of identification was needed. 

In 1968, two independent studies 36 * 37 of 
nonMnear minimization methods (output error 
methods) for obtaining aerodynamic coefficients 
were published, one describing the maximum likeli- 
hood estimator 36 * 3 ^ (with a Gauss-Newton tech- 
nique) to obtain a complete set of aerodynamic 
coefficients from flight data and the other 

describing a quasi 1 i neari zati on technique 37 * 39 to 
estimate some coefficients of an aircraft. One 
reason for the early success of these two methods 
is that previous research had furnished a well- 
defined model that adequately described the 
resulting motion of the vehicle. These two early 
results of aircraft identification by nonlinear 
minimization renewed interest in analysis of 
flight data. There was a later modification to 
these techniques to include a priori informa- 
tion.^ The minimization of this modified cost 
functional does not result in a maximum likelihood 
estimator, because it is based on the joint proba- 
bility distribution rather than the conditional 
probability. Other successful computer programs 
have been reported ,41-44 Extensive experience at 
many instal 1 ations 26 * 26 * 46-66 had been obtained 
using the maximum likelihood estimator technique 
on dynamic flight responses. 

Another approach, similar to these output 
error methods, was the application of the Kalman 
filter 69 to estimate the aerodynamic coefficients. 
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Some of the early results obtained by the Kalman 
filter technique were unsatisfactory; that is, the 
estimates of both the states and the parameters 
were biased and did not always converge to reason- 
able results. Improved results were obtained by 
adding the derivative of the state. 60 A weakness 
0 f the Kalman filter method is its dependence on 
the covariance matrix obtained from the filter. 60 
However, a technique was developed for obtaining 
estimates of the covariance matrix with a subop- 
timai Kalman filter. 6 9 A successful application 
of the Kalman filter to provide the state esti- 
mates used for the estimation of stability and 
control derivatives and performance parameters 
was subsequently described . 61 » 6 ^ 

Nondetermi ni St i c Analysis . As previously 
men t Toned” two classes of techniques were offered 
for the estimation of systems with measurement 
and state noise: the Kalnan filter (or more 

( jp nf » r a 1 |y 9 the extended Kalnan filter) tech- 
nique 18 * 1 9,^,60,63,64 and the maximum likeli- 
hood technique. 21 -23,65,66 The maximum likeli- 
hood estimator for the nondetermini Stic case is 
usually referred to as the filter error method. 

The general application of the extended Kalman 
filter was discussed in Refs. IB and 19. The 
extended Kalman filter for the discrete-time case 
was applied to simulated aircraft data with a 
state noise input.. 60 A similar application 63 to 
aircraft flight response data gave inconclusive 
results because the state noise input was small 
and the system was nonlinear. Somewhat better 
results were obained with an application of a 
greatly simplified extended Kalman filter 

t echo i que . 64 

The maximum likelihood estimator was applied 
to response data of an aircraft flying in atmos- 
pheric turbulence; 23 the resulting coefficients 
were in agreement with results obtained for the 
same aircraft flying in smooth air, that is, 
without state noise. 

Most of the results presented in this paper 

are based on an output error method program; 67 
the 1 1 i f f -Ha me code of this program is capable 
of using the Maine-Iliff formulation 68 (which 
can account for effects of state noise), although 
this feature is not used for the examples in 
this paper. 

basic Uses of Flight-Determined Coefficients 

The extraction of unknown aerodynamic coef- 
ficients or stability and control derivatives 
from flight data has been of interest for many 
years. The coefficients are used to provide 
final verification of the predicted full-scale 
design and to assist in the flight testing and 
verification of overall aircraft system perform- 
ance. 1 * 46 After the analysis of the flight 
tost data, the aircraft coefficients can be com- 
pared with calculated coefficients, estimates from 
computational fluid dynamics, and wind tunnel pre- 
dictions, and these comparisons can be used to 
update prediction methods for the improvement of 
future aircraft designs. 1 * 46 Once an aircraft is 


built, the coefficients play an important role in 

the expansion of the flight test envelope. 1 * 47 
As estimates of the derivatives become available, 
they are used to upgrade fixed-based simulators to 
assist in flight planning and aircraft control 
system modi f ication. T > 48 In addition, the flight- 
determined coefficients can be used to establish 
compliance with the desired design specifications. 
Flight-determined coefficients are also used to 
establish the accuracy of airborne simul ati ons 49 
and to identify aircraft parameters for adaptive 

control . 23 

Defin ition of Estim ation Technique 

The parameter estimation problem can be 
defined quite simply in general terms. The sys- 
tem under investigation is assumed to be modeled 
by a set of dynamic equations containing unknown 
parameters. To determine the values of the 
unknown parameters, the system is excited by a 
suitable input, and the input and actual system 
response are measured. The values of the unknown 
parameters are then inferred based on the require- 
ment that the model response to the given input 
match the actual system response. When formulated 
in this manner, the unknown parameters can be 
identified easily by many methods; however, com- 
plicating factors arise when application to a 
real system is considered. 

The first complication is the impossibility of 
obtaining perfect measurements of the response of 
any real system. The inevitable sensor errors are 
usually included as additive measurement noise in 
the dynamic model, and the theoretical nature of 
the problem then changes drastically. It becomes 
impossible to identify exactly the values of the 
unknown parameters; instead, the values must be 
estimated by some statistical criterion. The 
theory of estimation in the presence of measure- 
ment noise is relatively st raight forward for a 
system with discrete time observations, requiring 
only basic probability. 

The second complication of real systems is the 
presence of state noise. State noise is random 
excitation of the system from unmeasured sources, 
the standard example for the aircraft stability 
and control problem being atmospheric turbulence. 

If state noise is present and measurement noise is 
neglected, the analysis results in the regression 

al gori thm. 6 

When both state and measurement noise are con- 
sidered, 68 the problem is more complex than in the 
cases that have only state noise or only measure- 
ment noise. 

The final complication for real systems is 
modeling. It has been assumed throughout this 
discussion that for some value (called the best 
value) of the unknown parameter vector, the system 
is correctly described by the dynamic model. 
Physical systems are seldom described exactly by 
simple dynamic models, so the question of modeling 
error arises. No comprehensive theory of modeling 
error is available. The most common approach is 
to ignore it; any modeling error is simply treated 
as state noise or measurement noise, or both, in 
spite of the fact that the modeling error may be 
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m. 


deterministic rather than random. The assumed 

thp S rnnr dt h S J’ CS c * n then be ad j usted to include 
the contribution of the modeling error. This pro- 
cedure is not rigorously justifiable, but combined 
with a carefully chosen model, it is probably the 
best approach available. * 

It is possible to make a more precise, mathe- 
matically probabil istic statement of the parameter 
estimation problem. The first step is to^eHne 
e general system model (aircraft equations of 

~ • *’"> <«" <* »-<«" <» th.’.oJSo;,. 

discrete form as 

*(t()> - xy m 


^t) - f[x(t), u(t), 5] ♦ F(On(t) 


x*(t) * Ax(t) ♦ Bu(t) ♦ Fn(t) (6) 

2 (ti) = Cx ( t-j ) + Du(ti) + Gni (7) 

nn^taf’ B, C ’ 3nd ° are s >' stem matrices. For the 
no-state-noise case, the ?-(tj) term 0 f Fa m' 
can be approximated by * 

J C<to) = xo(0 (8) 

X * (t 1*l> - ♦xe(tf) ♦ *[u(ti) ♦ u(t i+1 )]/2 (9) 

^^(fi) - Cx^ft-j ) + Du(t-j) (10) 

the r rr»n!J rdnSiti0n matrix * and the integral of 
the transition matrix, *, are given by 


z(t i ) = 9 tx(tj ), u(t, ), 5] + G(?)hi (3) 

where x is the state vector, z i s the observation 
vector, f and g are system state and observation 
functions, u is the known control input vector 
C is the vector of unknown parameters, n is the 
state noise vector, n is the measurement noise 
*^ tor ’ F and G are system matrices, t is time, 

• denotes derivative with respect to time. 

The state noise vector is assumed to be zero-mean 

no x! Gau f s,an and stationary, and the measurement 
vector is assumed to be a sequence of inde- 
pendent Gaussian random variables with zero mean 
and identity covariance. For each possible Tstl 

Thf L’! the u,1 * n0wn parameters, a probability 
that the aircraft response time histories attain 
values near the observed values can then be 
defined. The maximum likelihood estimates are 
e ined as those that maximize this probability. 
Maximum likelihood estimation has many desirable 
statistical characteristics; for example, it 
yielos asymptotically unbiased, consistent, and 
efficient, estimates. 38 

,,, ,‘ f h th ! re no state no ’ se > then the maximum 
likelihood estimator m’nimizes the cost function 


♦ = exp[A(t i+1 - ti)] (ll a) 

f ti+1 

* • exp(Ai) dx B (ll b) 

Jz i 

basecTon thfn is J m P ortant . the estimator 

based on the nonlinear form of Eqs. (1) to (3) is 

intractable, and ad hoc techniques are required . 69 
,^,I 0 ^l n :. m1 f e tbe cost fu " c ”°n JU), we can 


app y the Newton-Raphson algorithm or some other 
minimization technique), which chooses successive 
estimates of the vector of unknown coefficients 


n k den ^ ing estimate) . if L is the iteration 
number, then the L + 1 estimate nf t ic „k. 
from the L estimate as eSt,mate of « « obtaine, 


1 " 

■MO = - X Lz(tj) - Z f (ti )]*(GG*)’ 1 


«ui • €L - Cv|j{i L )3-i[v*jfi L , J {l2] 

."wS; :n\n ss «* *»» <■ 


V (5) * * ■ * e (ti)]*(G 6 *) - 1 tv 5 2 c ( tl ,] 


f?(ti) - 2 e (tj )] + J N 1 n I (GG- 


7 l JU) = .itVc (ti)] * (GG *)" 1 tVr( t i)] 

1 * 1 s s 


where G(i* is the measurement noise covariance 
matrix, ?$( t, ) is the predicted response estimate 
of z at ti for a given value of the unknown- 
parameter vector £ (with ~ denoting predicted 
estimate), N is the number of time points, and 
* denotes transpose. The cost function is a func- 

r^n ? h T d,Fference between the measured and 
computed time histories. 

If Eqs. {?) and (3) are linearized (as is the 
and C ° ntr01 der,yat, - S in 

x ( t 0) = XQ 151 


' * 2 £( t i)l*(GG*)‘ 1 [ 7 ^(t i )] 


The Gauss-Newton approximation to the 
dient is 


^J(C) 2 2 !CV z (ti ) 3 *(GG*)- 1 [V z (ti 
1=1 s 4 £ € 


second gra- 


)] (14b) 


The Gauss-Newton approximation is 


much easier than the Newton-Raphson^ethod 1 because 
the second gradient of the innovation never needs 
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to he calculated. In addition, it can have the 
advantage of speeding the convergence of the 
algorithm, as is discussed in Ref. 6 . 


Figure I illustrates the maximum likelihood 
estimation concept. The measured response is com- 
pared with the estimated response, and the dif- 
ference between these responses is called the 
response error. The cost functions of Eqs. (4) 
and (11) include this response error. The minimi- 
zation algorithm is used to find the coefficient 
values that minimize the cost function. Each 
iteration of this algorithm provides a new esti- 
mate of the unknown coefficients on the basis of 
the response error. These new estimates of the 
coefficients are then used to update values of the 
coefficients of the mathematical model, providing 
a new estimated response and therefore a new 
response error. The updating of the mathematical 
model continues iteratively until a convergence 
criterion is satisfied. The estimates resulting 
from this procedure are the maximum likelihood 
estimates. 


The maximum likelihood estimator also pro- 
vides a measure of the reliability of each esti- 
mate based on the information obtained from each 
dynamic maneuver. This measure of the reliabil- 
ity, analogous to the standard deviation, is 
called the Cramer-Rao bound 5 * 24 or the uncer- 
tainty level. The Cramer-Rao bound as computed 
hy current programs should generally be used as a 
measure of relative accuracy rather than absolute 
accuracy. The bound is obtained from the approx- 
imation to the information matrix, H, which Is 
based on Eg. (14b); the actual information matrix 
is defined when evaluated at the correct values 
(not maximum likelihood estimates) of all the 
coefficients. The bound for each unknown is the 
square root of the corresponding diagonal ele- 
ment of H* 1 ; that is, for the ith unknown, the 
Cramer-Rao bound is\rH *{i,i). 


The formulation and the minimization algorithm 
previously discussed (Eqs. (4) to (14)) are imple- 
mented with the Iliff-Maine code (MMLE3 maximum 
likelihood estimation program). The program and 
computational algorithms are described fully in 
Ref. 67. All the computations shown and described 
in the remainder of this paper use the algorithms 
exactly as described in Ref. 67. 




with the wing in the forward position, and the 
T-37 aircraft at low speed. 

Derivation of an equation describing this 
motion is straightforward. Figure 2 illustrates 
an aircraft with the x axis perpendicul ar to the 
plane of the figure (positive forward on the 
aircraft). The rolling moment L‘, roll rate p, 
and aileron deflection 5 are positive as shown. 
For this example, the only state is p, and the 
only control is 5. The result of summing 
moments is 


I x p = L*(p, 6 ) (15) 

where I x is the rolling moment about the 
subscripted (x) axis. The first-order Taylor 
expansion then becomes 

p = Ik dp + d« (16) 

3p 3{ 

assuming small perturbations and using the nota- 
tion 


p = Lpp + Ls6 (17) 

where 


L = L7I X 

and the subscripts p and 5 denote partial deriva- 
tive with respect to the subscripted variable. 

Equation (17) is a simple aircraft equation 
where the forcing function is provided by the 
aileron and the damping by the damping-in-roll 
term Lp. In subsequent sections we examine in 

detail the parameter estimation problem where 
Eq. (17) describes the system. For this single- 
d eg ree-of- freedom problem, the maximum likelihood 
estimator is used to estimate L p or I 5 , or both, 
for a given simulated time history. 

We will assume that the system has measurement 
noise but no state noise; therefore, we can use 
Eqs. (1) to (3). Equation (4) then gives the cost 
function for maximum likelihood estimation. The 

weighting (GG*)~* is unimportant for this problem, 
so let GG* * 1. For our example. 


Simple Simulated Example 

For the discussion that follows, some knowl- 
edge of differential equations is assumed. A 
full derivation and a discussion of the aircraft 
equations of motion are given in Ref. 6 . 

The basic concepts involved in a parameter 
estimation problem will be illustrated by a simple 
simulated example representative of a realistic 
problem: an aircraft that exhibits pure rolling 

motion from an aileron input. This example, 
although simplified, typifies the motion exhibited 
by many Aircraft in particular flight regimes, 
such as the F - 14 aircraft flying at high dynamic 
pressure, the F-lll aircraft at moderate speed 


zi = xi 

Therefore, Eq. (4) becomes 

1 N ? 

J (Lp « L$) * — X [pi - Pi(Ep, Is)] (19) 

2 i=l 

where p«f is the value of the simulated measured 
response p at time ti and pi(Lp, L 5 ) is the esti- 
mated time history of p at time ti for L p = L p 
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and L$ = Lg. Throughout the rest of this paper, 
where simulated data (not experimental flight 
data) are used, the simulated measured time 
history refers to p^ , and the estimated computed 
time history, which varies with each iteration, 
is pi(Lp, 1 $). The estimated time history is 
a function of the current estimates of Lp and Lg, 

but the simulated measured time history, , 
is not. 

The most st rai ght forward method of obtaining 
Pi is with Eqs . { 8 ) and (9), Using the previously 
stated notation. 


Pi +1 = *Pi ♦ *(fii + Sj+i)/2 (20) 


where 

♦ = exp(lpA) 


( 21 ) 


* 



exp(Lpt) dT L ,5 


L{[1 - exp(LpA)] 
L P 


( 22 ) 


and A is the length of the sample interval, 
t i + 1 ~ t j • 


The maximum likelihood estimate is obtained by 
minimizing the cost function (Eq. (19)), which is 
done by applying the Gauss-Newton method. Equa- 
tion ( 12 ) is used to determine successive vof the 
estimates of the unknowns during the minimization. 

For this simple problem, £ = [Lp Lg]*, and 

successive values of L p and Lg are determined by 

updating Eq . (12). The first and second gradients 
of Eq. (12) are defined by Eqs. (13) and (14b). 

We now can write the entire procedure for 
obtaining the maximum likelihood estimates for 
this simple example. To start the algorithm, 
initial estimates of Lp and Lg are needed 

to define the value £g. Using Eq. (12), £ 1 and 

subsequently £^ are defined by using the first 

and second gradients of J(L p , Lg ) from Eq. (19). 

The gradients for this particular example, from 
Eqs. (13) and (14b), are 


N 

v (c l > = - Z (Pi - Pi) Vi (23) 

i =1 * 


O „ N 

V<e L ) s 2 (’ c Pi)*V c Pi (24) 

Computat ional Details of Minimization 


In the previous section we specified the 
equations for a simple example and described the 
procedure for obtaining estimates of the unknowns 
from a dynamic, maneuver. In this section we give 


the computational details for obtaining the esti- 
mates. Some of the basic concepts of parameter 
estimation are best shown with simulated measured 
data, where the best (correct, in this simulated 
case) answers are known. Therefore, in this sec- 
tion we study two examples involving simulated 
time histories. The first example is based on 
data that have no measurement noise, which results 
in estimates that are the same as the correct 
values. The second example contains significant 
measurement noise; consequent ly , the estimates are 
not the same as the correct values. 

For this simulated example, 10 points (time 
samples) are used. The simulated measured data, 
which we refer to as the measured data, are based 
on Eq. (17). We use the same correct values 
L p ~ -0.2500 and Lg = 10.0 for both examples. 

In addition, the same input 6 is used for both 
examples, the sample interval A = 0.2 sec, and 
the initial conditions are zero. Tables of all 
the significant intermediate values of the calcu- 
lations are given in Ref. 6 .^ In both examples, 
the initial values defining £ are L D = -0.5 and 
L <5 = 15.0. u 

Example With No Measurement Noise . The simu- 
lated measured time history of aileron deflection 
for the case with no measurement noise (no-noise 
case) is shown in Fig. 3. The aileron input 
starts at zero, goes to a fixed value, and then 
returns to zero. The resulting simulated measured 
roll rate time history is also shown. 

Table 1 gives the values for L p , Lg, and J for 

each iteration, along with the values of <t> and tp 
needed for calculating p-j . In three iterations 
the algorithm converges to the correct values to 
four significant digits for both L p and Lg. 

Figure 4 shows the match between the simulated 
measured data and the estimated data for each of 
the first three iterations. The match is very 
close after two iterations and is nearly exact 
after three. 

Although the algorithm converges to four-digit 
accuracy in L p and Lg, the value of the cost func- 
tion J continues to decrease rapidly between 
iterations 3 and 4. This is a consequence of 
using the maximum likelihood estimator on data 
having no measurement noise. Theoretically, with 
infinite accuracy the value of J at the minimum 
should be zero. However, with finite accuracy the 
value of J becomes small but never reaches zero. 
This value is a function of the number of signifi- 
cant digits. For the 13-digit accuracy used here, 
the cost eventually decreases to approximately 
O.cS x 10-23. 

Example With Measurement Noise . The simulated 
measured data used in the case with measurement 
noise (noisy case) are the same as those used in 
the previous section, except that pseudorandom 
Gaussian noise is added to the roll rate (Fig. 5). 
The si gnal -to-noi se ratio is quite low in this 
example (compare Figs. 3 and 5). The values of 
L p , Lg, <p, ip, and J for each iteration are given 
in Table 2. The algorithm converges in four 
iterations. The behavior of the coefficients as 
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they approach convergence is much like that in the 
no-noise case. The most notable result of this 
case is that the converged values of Lp and L$ are 
somewhat different from the correct values. The 
natch between the simulated measured and estimated 
time histories is shown in Fig. 6 for each itera- 
tion. No change in thr match is apparent for 
iterations 2 and 3. The match is very good con- 
sidering the amount of measurement noise. 

In Fig. 7, the tine history estimated using 
the no-noise estimates of L p and L$ is compared 
with that using the noisy estimates of L p and L$. 
Because the algorithm converged to values somewhat 
different from the correct values, the two esti- 
mated time histories for their respective values 
are similar but not identical. 

The accuracy of the converged elements can be 
assessed by looking at the Cramer-Rao inequal- 
ity? 4 , 67 discussed previously. The Cramer-Rao 
bound can be obtained from an approximation to 
the information matrix H, where 

( N |-1 

»* = ?J min| X [V c 2 c (ti)l*{GG*)' 1 V 5 2 e {t i )| /(N-l) 

The Craner-Rao bounds for L p and L$ are the square 

roots of the diagonal elements of the H"* matrix, 

or Vh ^ (1,1) and Vh _1 (2,2), respectively. The 
fram^r-Rao bounds are 0.1593 and 1.116 for Lp and 

1$, respectively. The differences between Lp and 
L p and between L 5 and are less than their 
respective bounds. 

Cost ^u nc t 10ns 

In the previous section we obtained the maxi- 
num likelihood estimates for simulated time 
histories by minimizing the values of the cost 
functions. To fully understand what occurs in 
this minimization, we must study in more detail 
t.he form of the cost functions and some of their 
more important charact.eri st i cs . In this section, 
the cost function for the no-noise case is dis- 
cussed briefly, The cost function for the noisy 
case is then discussed in more detail. The same 
two t»ne histories studied in the previous section 
are examined here. The noisy case is more inter- 
esting because it has a meaningful Cram§r-Rao 
hound and is more representative of aircraft 
flight, data. 

It is important to remember that in this paper 
everything related to cost functions (Eq. (19)) is 
based on simulated tine histories that are defined 
hy Eq. (17). For every measured time history we 
night choose (simulated or flight data), a com- 
plete cost function is defined. For the case of 
n variables, the cost function defines a hyper- 
surface of n + 1 dimensions. We could avoid 
bothering with the minimization algorithm if we 
could construct this surface and look for the 
minimum, but this is not a reasonable approach, 
because the number of variables is generally 
greater than two. Therefore, the cost function 
can be described mathematically but not pictured 
graph! ca 1 ly . 


rfcV • ** 


One-Dimensional Case . To illustrate the many 
aspects of cost functions, it is easiest to look 
first at cost functions having one variable. In 
an earlier section, the cost function of L p and L$ 
was minimized. That cost function is most inter- 
esting in the L p direction. Therefore, the one- 
variable cost function studied here is J(L p ), 
with the correct value of L 6 = 10. U. Figure 8 
shows the cost function plotted as a function of 
L p for the no-noise case. As expected for this 
case, the minimum cost is zero and occurs at the 
correct value of L p = -0.250U. It is apparent 
that the cost increases much more slowly for a 
more negative L p than for a positive L p . In fact, 

the slope of the curve tends to become less nega- 
tive where L p < -1.0. Physically this makes sense 

because the more negative values of L p represent 
cases of high damping and the positive L p repre- 
sents an unstable system. Therefore, the pj for 
positive Lp becomes increasingly different from 
the measured time history for small positive 
increments in L p . For very large damping (very 
negative L p ), the system would show essentially no 
response. Therefore, further large increases in 
damping result in relatively small changes in the 
value of J(L p ) . 

In Fig. 9, the cost function based on the 
noisy case time history is plotted as a function 
of L p . The correct L p value (-0.2500) and the 

L p value (-0.3218) at the minimum of the cost 
(3.335) are both indicated on the figure. The 
general shape of the cost function in Fig. 9 is 
similar to that shown in Fig. 8, Figure 10 com- 
pares the cost functions based on the noisy and 
no-noise cases. The comments relating to the cost 
function based on the no-noise case also apply to 
the cost function based on the noisy case. Figure 
10 shows clearly that the two cost functions are 
shaped similarly but shifted in both the L p and J 

directions. Only a small difference in the value 
of the cost would be expected far from the minimum 
because the "estimated" time history is so far 
from the simulated measured time history that it 
becomes irrel event as to whether the simulated 
measured time history has noise added. Therefore, 
for large values of cost, the difference in the 
two cost functions should be small compared with 
the total cost. 

Figure 11 shows the gradient of J(L p ) plotted 
as a function 0. L p for the noisy case. Finding 
the zero of this function (or equivalently, the 
minimum of the cost function) using the Gauss- 
Newton method was discussed previously. The 
gradient is zero at L p = -0.3218, which corre- 
sponds to the value of the minimum of J(L p ). 

The usefulness of the Cramer-Rao bound was 
discussed in the Example With Measurement Noise 
section. It is useful to digress briefly to 
discuss some of the ramifications of the Cramer- 
Rao bound for the one-dimensional case. The 
Cramer-Rao bound has meaning only for the noisy 
case. In the noisy example, the estimate of L p 

is -0.3218, and the Cramer-Rao bound is 0.0579. 
The calculation of the Cramer-Rao bound was 
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Two-Dimensional Case. In this section, the 
rost”” function depi^Ton both L p and Lj is 
studied. The no-noise case is examined first, 
followed by the noisy case. 

Even though the cost function of 

only two unknowns, it. is much more d fficult * 
visualize than is the one-dimensional case. The 


cost function over reasonable ranges of Lp and L« 

• c-in i2 The minimum must lie in the 
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nearly quadratic there, so fairly rapid conv - 
gence in this region would be expected The L p 
axis becomes an asymptote for cost as L« approach 
zero. The cost is constant for Ls • 0 because 
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Lp, resulting in constant cost. 
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U 9 = - 0.2500 and U = 10.0. This is also evident 
Kw innkina at the cost function surface shown in 
Fig! 14. The surface has its minimun. at the cor- 
rect value. As expected, the value of the cost 
function at the minimum is zero. 

As in the one-dimensional case, the primary 

the minimum was at a mgner tuai a 
Live value of L p . In the two-dimensional case, 
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follow the path of the minimization example used 
before. The first iteration (L = 1) brought the 
values of Lp and L$ very close to the values at 
the minimum, and the second essentially arrived at 
the minimum (viewed at this scale). One of the 
reasons the convergence is so rapid in this region 
is that the isoclines are nearly elliptical, dem- 
onstrating that the cost function is very nearly 
quadratic in this region. If we had started the 
(iauss -Newton algorithm at a point where the iso- 
clines are much less elliptical (as in some of the 
border regions in Fig. lb), the convergence would 
have progressed more slowly initially, but it 
would have progressed at much the sane rate as it 
entered the nearly quadratic region of the cost 
tunc.t ion. 

Before concluding our examination of the two- 
dimensional case, we shall examine the Cram6r-Rao 
bound. Figure 18 shows the uncertainty ellipsoid, 
which is based on the Cram6r-Rao bound. (The 
relationships between toe Cramer-Rao bound and the 
uncertainty ellipsoid are discussed in Ref. 69.) 

The uncertainty ellipsoid almost encloses the 
correct values of Lp and Lg. The Cram£r-Rao bound 

for Lp and L$ can be determined from the projec- 
tion of the uncertainty ellipsoid onto the Lp and 
Lfi axes and then compared with the values calcu- 
lated for the noisy case, which were 0.1593 and 
1.116 for Lp and Lg, respectively. This projec- 
tion is analogous to the case for n unknowns, but 
in that case the projection would be the n ♦ 1 
hyperellipsoid's projection onto a hypersurface. 

Estimation Using Flight Data 

We have examined the basic mechanics of 
obtaining maximum likelihood estimates from sim- 
ulated examples with one or two unknown param- 
eters. To make the transition from theory to 
practical application, we present results obtained 
from analysis of actual flight data and discuss 
how the aircraft parameter estimation results are 
used to solve real problems. In this case we 
illustrate the necessity of obtaining estimates of 
the aircraft coefficients of the differential 
equations of motion (the stability and control 
derivatives) to solve important and related prob- 
lems encountered in flight. However, the aircraft 
stability and control example is only one of sev- 
eral applications of parameter estimation tech- 
niques; useful results can be obtained in many 
applications where the phenomenology is well 
understood. For the computat i onal 1 y difficult 
situation usually encountered with actual flight 
data, we obtain the maximum likelihood estimates 

with the Iliff-flaine code (MMLE~ program ). ^ 

Before studying the specific examples, a brief 
historical review of some other uses of the esti- 
mates is presented. 

In the past, the primary reason for estimating 
stability and control derivatives from flight 
tests was to make comparisons with wind tunnel 
estimates. As aircraft became more complex and as 
flight, envelopes were expanded to include flight 
regimes that were not well understood, new re- 
quirements of the derivative estimates evolved, 
ior many years, the f 1 i ght-determi ned derivatives 
werp used in simulations to aid in flight planning 


and in pilot training. The simulations were par- 
ticularly important in research flight test pro- 
grams in which an expansion of the envelope into 
new flight regimes was required. As more was 
learned about these new flight regimes, the com- 
plexity of the aircraft, and particularly their 
sophisticated flight control systems, increased. 

The design and refinement of the control system 
for these complex aircraft required higher fidel- 
ity simulations. As a consequence, a more com- 
plete knowledge of the flight-determined stability 
and control derivatives was necessary. Almost 
all current high-performance aircraft have very 
complex control systems to compensate for their 
deficiencies in basic aerodynamic characteristics. 
Consequently, most flight test programs for these 
aircraft require a complete flight-determined 
set of stability and control derivatives, and 
parameter estimation techniques for estimating 
stability and control derivatives from flight 
data have become more sophisticated. 

At the Dryden Flight Research Facility of 
NASA's Ames Research Center ( Ames- Dryden ) , ana- 
lysts have been involved in the estimation of sta- 
bility and control derivatives with maximum like- 
lihood estimators since 1966 and have success- 
fully applied maximum likelihood estimators to 
nearly 50 different aircraft configurations. Some 
of the experience gained through these applica- 
tions Is included in che bibliography of Ref. 9. 
Recent Ames-Dryden applications have concentrated 
on estimating stability and control derivatives 
to assist in designing or refining control sys- 
tems. Three such applications (to be discussed in 
detail) are the F-14, highly maneuverable aircraft 
technology (Hi MAT), and space shuttle programs. 

All three of these programs have made extensive 
use of high-fidelity, pil ot-in-the-loop simula- 
tions, which are implemented using the best wind 
tunnel data available. Portions of these flight 
test programs were defined to obtain data for 
refining simulator models. 

The chosen method of enhancing the simulator 
model depends on the aircraft involved in the 
flight test program. The F-14 aircraft flew sev- 
eral flights specifically for defining the sta- 
bility and control derivatives over a large angle- 
of-attack range because the necessary control 
refinement related to the hi gh-angl e-of-attack 
regime. The HiMAT vehicle flew several flights 
with a positive static margin (stable open-loop 
system) so that derivatives could be obtained to 
design a control system for flight at a negative 
static margin (unstable open-loop system). The 
space shuttle entered from space on the most con- 
servative trajectory to allow assessment of its 
characteristics before an envelope exoansi on 
was begun. 

Once the flight data are obtained and ana- 
lyzed, the simulator is updated to assist in con- 
trol system design and further flight planning. 
Where flight results agree with wind tunnel pre- 
dictions, confidence in the simulation grows, and 
envelope expansion proceeds more efficiently. 

The coefficients evaluated in this section 
are contained in the aircraft equations of motion, 
which are derived and discussed in detail in 
Ref. 6. 
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F -14 Ai rcraft 


The F-14 aircraft is a twin-engine, high- 
performance fighter with variable wing sweep 
(Fig. 19). The Ames-Dryden F —1 4 program was 
intended to improve the handling qualities of the 
airplane at high angles of attack by incorporating 
several conrr,.] system techniques . 7U, 71 j he f i rst 
part of the program was dedicated to obtaining 
flight -determined stability and control deriva- 
tives for the subsonic envelope of the F-14 air- 
craft, the complete trimmed angl e-of-attack range 
for Mach number M * 0.9. 

In many instances the flight data agreed 
with the wind tunnel predictions; Fig. 2i) (from 
Ref. 70) shows the comparison of C n(3 (C n being the 

coefficient of yawing moment) as a function of 
angle of attack a from flight and wind tunnel 
estimates. (Througout this and following dis- 
cussions, a subscript to the coefficient denotes 
partial derivative with respect to the subscripted 
variable.) The symbols denote the estimate, and 
the vertical bar designates the uncertainty level 
(Cramer-Kao bound). The agreement is good, 
although there is some disagreement at a >'25°; 
nevertheless, the same trends are seen for both 
flight and wind tunnel data. 

Figure 21 shows the fl ight -determined C 

*P 

(C t being the coefficient of rolling moment) 

as a function of a for M < 0.55 and for M = 0.9. 
There was some uncertainty in the accuracy of the 
wind tunnel predictions of C ^ because the wind 

tunnel model configuration was different from the 
flight conf i guration. The implementation of 
c fcp at M = in the simulation produced a pre- 
viously unsimulated wing rock characteristic that 
had been observed in flight. The wing rock had 
been a troublesome characteristic, and its simula- 
tion was important in improving handling qualities 
through control system modifications. Figure 22 
shows the flight-determined values of C 0 as a 
, 

function of a compared with the results of two 
different sets of wind tunnel results. There had 
been some concern about the disagreement between 
the two sets of wind tunnel results before flight. 

At low angles of attack, the three sets of esti- 
mates are in fair agreement; however, at a > 15°, 
the flight data lie between the two sets of wind’ 
tunnel data. 

A last example from the F-14 aircraft shows 
how the wind tunnel and flight estimates inter- 
play to improve a simulation. After the lateral- 
directional derivatives were incorporated in the 
simulation, the resulting simulated lateral- 
directional motions from a longitudinal -stick snap 
maneuver were found to be inconsistent with the 
flight response. Since the F-14 program was pri- 
marily a lateral-directional investigation, the 
longitudinal derivatives in the simulation had not 
been updated with the flight-determined values. 

When the flight-determined longitudinal deriva- 
tives were included in the simulation, the stick 
snap response agreed more closely with the flight 
response. In tracking down the inconsistency, a 
large discrepancy was discovered between the wind 


tunnel and flight-determined values of C m (C m 

being the coefficient of pitching momont). This 
is shown in Fig. 23, where fl ight -determined 
C ma is compared with the wind tunnel estimates 

of c m a for the untri [timed and trimmed conditions. 
Further investigation showed that the untrinmed 
values of C m(j had been put in the simulation and 

that the predicted trimmed values of C m were in 
excellent agreement with flight estimates. 


f 1 ight data , i n add 
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a large part of the envelope. Since the 
F-14 hi gh-ang I e-of-attack program also needed to 
examine responses of a highly transient nature, 
more tedious and time-consuming fine tuning of the 
simulation was required for flight at other than 
near the trimmed conditions. With the resulting 
simulation, the proposed control system techniques 
were further refined; the result was a more effi- 
c lent demonstration in flight. 

This exemplifies the value of flight test 
parameter estimation in improving the handling 
qualities of an aircraft through control system 
improvements. 


Hi MAT Vehi cl e 


The HiMAT vehicle is a remotely piloted 
research vehicle with advanced close-coupled 
canards, wing-type winglets, and provisions for 
variable leading-edye camber. It is made of 
advanced composite materials to allow for aero- 
elastic tailoring and to minimize weight. It 
was flown in an unstable configuration because 
the wing deformation then resulted in a desirable 
camber shape at high load factor and because the 
trim drag was reduced. 








. X a- UCb I yiltfU 

to fly With a sustained 8-g turn capability at 
Mach 0.9 at 25,000 ft altitude and to demonstrate 
flight supersonically to Mach 1.4. To attain the 
Mach 0.9 condition, it was predicted that the 
vehicle must be flown in an unstable configurator 
(10-percent mean aerodynamic cord (MAC) negative 
static margin). The philosophy for testing the 
HiMAT vehicle was somewhat different from that for 
production aircraft: Flight-determined stability 
and control derivatives were to be relied on to ' 
keep the wind tunnel program to a minimum. The 
original simulation data base contained the wind 
tunnel data supplemented with some computed 
characteristics. 


me vemcie was flown in 


— stabl e conn g- 
uration to obtain stability and control deriva- 
tives with the control feedbacks set to zero. 
While these data were being gathered, a control 
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system suitable for unstable flight was being 
designed* based on wind tunnel tests. Then, with 
the fl ight-determi ned derivatives, the simulator 
was updated and the control system adjusted for 
this update so' that the unstable vehicle could be 
flown safely. Stability and control maneuvers 
were performed at subsonic and supersonic Mach 
numbers, at angles of attacx up to 10°, and at 
altitudes from 15,000 to 45,000 ft. A complete 
set of stability and control characteri sti cs was 
obtained for both the longitudinal and lateral- 
direciional degrees of freedom. 75 Because the 
values of the HiMAT derivatives are classified, 
the data are plotted on unlabeled vertical axes; 
nevertheless, an assessment of predicted and 
f 1 ight -determi ned derivatives can still be made. 
All the derivatives, predicted and flight deter- 
mined, are corrected to 0-percent MAC. For the 
flight conditions flown, there were no aeroelastic 
effects noted in the flight data. 


Figure 25 shows fl ight-determ ned directional 
dynamic stability as a function of Mach 

number at a = 4° compared with the rigid and flex- 
ible predictions. Flight estimates are about the 
same as predictions at M = 0.4 and 0.9, but they 
differ significantly in between. In Fig. 26, 


n Sdyn 


is plotted as a function of a at M - 6.9, 


showing that the vehicle is slightly unstable in 
the lateral -di rectional axes at the lower angles 
of attack. Considering that these data are 
plotted for 0-percent MAC , this instability would 
be considerably greater and over a wider angle-of- 
attack range if the center of gravity were moved 
significantly aft (aft movement of the center of 
gravity makes any vehicle less stable). The de- 
rivatives C y „ ( C y being the coefficient of side- 

D 


force) and C. agreed with predictions; however, 
Cn r was twice the predicted value, C np was of 
opposite sign, and C^ was a small fraction of 

predictions. The rolling moments due to aileron, 
C 9 , agreed fairly well with the rigid predic- 
* 6 DE 

tions; C nA was 25-percent less than predicted; 
o r 


both Cn 6 ^ and C n<5 ^ showed a positive increment 

over prediction. The derivative C { was about 

6 r 

twice the predicted value. Since there were so 
many large differences between the flight- 
determined derivatives and the minimal wind tunnel 
set, it was decided to completely reevaluate the 
1 ateral -di recti oral control laws designed for the 
unstable configuration using the flight data 
instead of the wind tunnel data, which we r e used 
in the original design. Some reasons for this can 
be seen in Fig. 27, in which the control deriva- 
tives Cn 6()E , Ct 6r# and Cp 6r are plotted func- 
tions of a at M = 0.9. These differences between 
flight and predicted values meant that the simu- 
lator had to be extensively revised. 


The HiMAT vehicle program was a technology 
demonstration program and therefore was required 
to demonstrate the technology only at specific 
design points. A technology demonstration is 
quite different from many programs, such as the 
F -14 program, because only certain steady-state 


requirements must be demonstrated. Therefore, all 
the points (or flight conditions) that needed to 
be flown were near steady-state points for which 
fl iyht -demonstrated derivatives already existed. 

To update the simulator, all the predicted data 
were disregarded, and only flight-determined sta- 
bility and control derivatives were used. The 
knowledge that the aircraft stability and control 
derivatives exhibited no significant aeroelastic 
effects permitted the reevaluation o f the unstable 
control system, and the design was simplified. 

The control laws designed for the unstable 
configuration were much more complex than the 
rate-feedback system used for gathering stability 
and control derivatives. The new control laws 
were modified by (1) adding a lateral acceleration 
ay feedback to improve closed-loop directional 
dynamic stability; (2) adding an interconnect 
between lateral stick and rudder to improve 
lateral control characteristics; (3) changing 
the various feedback gains' to improve damping 
characteristics; and (4) locking the aileron 
surface to eliminate adverse yaw and also to 
eliminate the possibility of a predicted surface- 
buzz problem at higher Mach numbers. This design 
of the lateral -di rectional control system was the 
result of an extensive study of possible control 
systems using both the simulator and the linear 
analysis techniques. When the new control system 
was designed, it was implemented on the HiMAT 
vehicle, and it was flown in a stable configura- 
tion. Control surface doublets were input, and 
the responses were compared with the simulator- 
derived responses. The comparison was excellent, 
giving confidence that the unstable vehicle could 
be tested. 

The benefits of flying the unstable vehicle 
were demonstrated in flight when a 0.4-g improve- 
ment in sustained-y capability was realized by 
changing the center-of-gravi ty location from the 
point of neutral stability to 5-percent MAC aft of 
the neutral point. When the unstable vehicle was 
flown with a 5-percent MAC negative static margin, 
a sustained turn of about 7.8 g was achieved. 

Based on these numbers, the HiMAT vehicle should 
be able to demonstrate a sustained 8.U-g turn 
capability with the 10-percent MAC negative static 
margin (unstable vehicle). 

In the case of the HiMAT vehicle, flight test 
parameter estimation became the sole method of 
defining the stability and control derivatives. 

A control system design for the unstable config- 
uration was defined from flight test results. 

The adequacy cf the design was demonstrated on 
the simulation updated with flight data. The 
resulting control system enabled the unstable 
vehicle to be flown. 

A recent investigation of determining the 
aerodynamic coefficients for the Siighly unstable 
X-29A vehicle is described in Ref. 69. This 
investigation sheds new light on parameter estima- 
tion of unstable systems, which has widespread 
application to systems other than those defined 
by stability and control derivatives. 

Space Shuttle urbiter 

The space shuttle orbiter is a large double- 
delta-winged vehicle designed to enter the atmos- 
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phere and land horizontally. The entry control 
system consists of 12 vertical reaction cont r ol 
system (RCS) jets (6 up-firing and 6 down-firing) 
and 8 horizontal RCS jets (4 left-firing and 
4 right-firing), 4 elevon surfaces, a body flap, 
and a split rudder surface (Fig. 28). The ver- 
tical jets and the elevons are used for both pitch 
d n,i roll control. The jets and elevons are used 
symmetrical 1 y for pitch control and asymmetric- 
ally for roll control. More information on the 
configuration and flight plan is given in Ref. 76. 

The F -1 4 and Hi MAT examples showed how param- 
eter estimation can be used in an incremental 
flight test program, that is, a progressive expan- 
sion of the flight envelope to obtain data in the 
'■iore certain areas first and in the more chal- 
lenging or hazardous ones later. However, the 
space shuttle program could not be approached in 
this manner, f or the vehicle had to demonstrate on 
the first flight that it could be flown safely 
over most of its envelope. Further complicating 
the program, this first flight included very 
hazardous flight regimes. The subsonic flight and 
landing characteristics had Deen demonstrated in 
the earlier approach and landing test program, 
but the hypersonic, peak heating, and transonic 
regions were largely unexplored for a vehicle of 
this typo. 

f.xtensive wind tunnel tests were performed, 
and those data were incorporated into high- 
fidelity simulations. No matter how carefully 
wind tunnel tests are performed, there are fre- 
quently di sc repanc i es between the predictions and 
the demonstrated flight character! st i »_s ; there- 
fore, uncertainties were defined for each stabil- 
ity and control derivative. These uncertainties 
(called variations in Rpt. 771 were based to a 
large extent on previously reported di screpancies 
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between predictions and flight.' 0 

In preparation for the first flight, a con- 
trol system was developed to provide satisfactory 
closed- loop vehicle characteristics fo r deriva- 
tives that fell between tne variations that had 
been previously defined. After flight data were 
obtained, the flight estimates of the stability 
and control derivatives were used to r educe the 
preflight variations. This reduction then allowed 
the control engineers to refine the control system 
and therefore to improve the shuttle handling 
qualities. ! n addition, tne flight-determined 
derivatives were used to determine if configura- 
tion placards (limitations on tne flight envelope) 
could be modified or removed. 


V = 24,300 ft/sec and at a dynamic pressure 

q - 12 1 b/ft 2. 

The actual STS -1 maneuver that occurred at 
this flight condition is shown in Fig. 30, which 
depicts a more hazardous maneuver than was pre- 
dicted. At this flight condition the excursions 
must be kept small. The flight maneuver resulted 
in twice the angle-of-sidesl ip 0 peaks predicted 
and in a somewhat higher roll rate than predicted. 
Also, there was more yaw-jet firing than was pre- 
dicted, and tne notion was more poorly damped than 
predicted. It is obvious from comparing the pre- 
dicted with the actual maneuver (Fig. 31) that the 
stability and control derivatives were signifi- 
cantly different than pred'cted. It is fortunate 
that the control system design philosophy dis- 
cussed previously had beei used for the shuttle. 
Although the flight maneuver resulted in excur- 
sions greater than planned, the control system 
did manage to damp out the oscillation in less 
than 1 min. With a less conservative design 
approach, the resulting entry maneuver could 
have been a good deal worse. 


To assess the problem with the first bank 
maneuver, the flight-determined stability and 
control derivatives were compared with the pre- 
dictions. Of all the derivatives obtained from 
STS-1, the two important ones that differed most 
from the predictions at the flight condition being 

discussed were C fl and the rollinq moment due to 
*6 

yaw jet firing, Lyj. Since the entry tends to 

monotonical ly decrease in Mach number, the deriva- 
tives can be best portrayed as functions of the 
guidance system "Mach number," which is V/1000. 
Figure 32 shows as a function of guidance Mach 

number, and Fig. 33 shows Lyj as a function of 

guidance Mach number. Only the estimates from 
STS-1 are shown in these figures. 


When inly the change in C. was entered into 

the simulation data base, the maneuver looked 
very much like the o r iginal prediction (Fig. 29); 
however, as expected, tne frequency of the oscil- 
lations changed to be more representative of the 
actual flight frequencies (Fig. 30). The effect 
on the simulation of changing only Lyj from the 


predictions is shown, with the flight response, in 
Fig. 34. These two tine histories are very close, 
considering that the other differences between 
the fl ight -determined and predicted derivatives 
have been ignored. 


Some of the stability and control results 
obtained from tne first three flights are con- 
tained in Refs. 79 and 80. One interesting 
example of wnere parameter estimation played an 
important role in the shuttle program occurred 
during the first energy -management bank maneuver 
on the first entry of the shuttle (STS-1). The 
response to the automated control inputs computed 
using the predicted stability and control deriva- 
tives is shown in Fiq. 29. It. should be noted 
that the control inputs shown here (and for all 
other simulation comparisons) are the closed- 
loop commands from the shuttle control laws. 

The maneuver was to be made at a velocity 


It is apparent that the primary problem with 
the initial bank maneuver was the poor prediction 
of Lyj. The control system software is very 

complex, and it cannot be changed and verified 
between shuttle missions; therefore, an interim 
approach was taken to keep this large excursion 
from occurring on future flights. The flight- 
determined derivatives were put into the simula- 
tion data base, and the shuttle pilots practiced 
performing the maneuver manually, trying to attain 
a smaller response within more desirable limits. 
The maneuver was performed manually on STS-2 to 
STS-4. Figure 36 shows the manually flown maneu- 
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vpr from STS-2. For this maneuver, roll rate, 
yaw rate, and sideslip angle were within the 
desired Units. The maneuver does not look like 
the original predicted response, because the de- 
rivatives and the input were different and the 
basic control system remained unchanged. Since 
the response variables were kept low and the inputs 
were slower and smaller, the flight responses on 
STS-2 to STS-4 did not show a tendency to oscil- 
late. The software was updated for STS-5, and 
thp resulting automated maneuver is essentially 
1 ndist i ngui shable from that shown in Fig. 35. 

This maneuver has been used on all subsequent 
shutt 1 e f 1 i qht s . 

The application of parameter estimation tech- 
mques to the highly complex space shuttle vehicle 
will continue, and the results of this application 
have and will si gni f icantly affect the control 
system design, placard modification, and flight 
procedures in general . 

Concluding Remarks 


^Schweppe, Fred C., Uncertain Dyna mic Systems, 
Prentice-Hall, Englewood Cliffs, New Jersey, 1973. 

5 Maine, Richard E., and II iff, Kenneth W., 

lajM #I? at10n ? f D *" amic Systems," AGARD-AG-300, 
1984. (Also published as NASA RP-1138, 1985.) 

^Maine, Richard E., and Iliff, Kenneth W., 
Application of Paraneter Estimation to Air- 
craft Stability and Control - The Output-Error 
Approach," NASA RP-1168, 1986. 

in B!n rt0n ’ F ,;^’ ' The Measu cement of the Damping 
Roll on a JN4h in Flight," NACA Rep. 167, 1923. 

^Norton, F.H., "A Study of Longitudinal 
Dynamic Stability in Flight," NACA Rep. 170, 

9 Iliff, Kenneth W., and Maine, Richard E., "A 
Par “ et ' r 


In this paper, the aircraft parameter estima- 
tion problem is used as an example of how param- 
eter estimation can be applied in many scientific 
and engineering fields to assess phenomenology 
from observations, and a brief survey of the 
literature is presented. The theory, a simple 
simulated example, and the application of experi- 
mental results to solvn real problems are given 
and explained. The maximum likelihood parameter 
estimation technique was used in the F -14 program 
to effect control system chanaes that improved 
handling qualities at high angles of attack. The 
same technique provided the primary source of 
information for control system refinement on the 
unstable Hi MAT vehicle. Space shuttle energy- 
management maneuvers have been redefined based on 
simulations using f 1 i ght-determined stability and 
control estimates. Moreover, parameter estima- 
tion techniques are being relied upon for future 
control system design, placard modification or 
removal, and flight procedures in general for 
the space shuttle. 

Thp explanation of parameter estimation tech- 
niques and the demonst rat i on of their highly suc- 
cessful application to the aircraft problems are 
intended to inform and to encourage scientists in 
other fields to consider these techniques for 
application to problems where a representat i ve 
model and high-quality data exist. 

References 


10 Gauss, Karl Friedrich, Theory of the Mot ion 
of the Heavenly Bodies Hoy i ng About the Sun in 
Ionic Sections (transl ated by fW l « Uenry 

?oa 1 7 S ’ H? 5 7 U Dover Publications, Inc., New York, 
lo4/. (Translated from Theoria Mot us . 1809.) 

^Douglas, J., "Theorems in the Inverse 
Problem in the Calculus of Variations,'* Proc. 
Natl. Acad. Sci., vol . 16, no. 3., 1940. 


^Gelfand, I.M. 
Determination of a 
Spectral Function," 
Nauk, USSR, 1951. 


, and Levitan, B.M., "On the 
Differential Equation From Its 
Izvestiya Ac., vol. 15, no. 4, 


l^Feldbaum, A. A., "Dual Control Theory," 
Automn. Remote Control, vol. 22, 1961, pp/l-12 
and 109-121. 


^Cuenod, M., and Sage, A., "Comparison of 
Some Methods Used for Process Identification *' 
Automatic a, vol. 4, 1968, pp. 235-269. 


lb Balakrishnan, A.V., and Peterka, V., 

" Identi fication in Automatic Control Systems " 
Automat l ca , vol. 5, 1969, pp. 817-829. 


l 6 Astrom, K.J., and Eykhoff, P. f "System 
Identi fi cation — A Survey," Automatica, vol. 7, 
no. 2, 1971, pp. 123-162. 


Mliff, Kenneth W., and Maine. Richard E., 
"NASA Dryden's Experience in Parameter Estimation 
and Its Use in Flight Test," AIAA-82-1 373, MAA 
AFM Conference, San Diego, California, Aug. 1982. 

2 Iliff, Kenneth W., and Maine, Richard E., 
"More Than You May Want to Know About Maximum* 
Likelihood Estimation," NASA TM-85905, 1985. 


^Eykhoff, Pieter, System Identification. 
Lonoon te ig % d State Est1mation » John Wiley * Sons, 

18 Astrom, Karl J., "Control Problems in 
Papermaking," Proc. IBM Scientific Computing 
Symposium on Control Theory and Applications, IBM 
Data Processing Division, White Plains, New York 
1966, pp. 135-167. 


^Zcdeh, Lot f i A., and Desoer , Charles A., 
^l^ -^stem Theo ry. McGraw-Hill Book Co., New 


^Kashyap, R.L., "A New Method of Recursive 
Estimation in Discrete Linear Svstems," IEEE 
Trans. Automatic Control, vol. AC-15, no 1 
Feb. 1970, pp. 18-24. ’ ’ 


14 


^’lidldkrishnan, A.V., "Stochastic System 
Identification Techniques," Stochastic 
Opt ini zat i on and Control , H.T7 Karreman, erf., 

John Wiley * Sons, London, 1968, 

^Balaknshnan, A.V., “Stochastic Differen- 
tial System l," Filtering and Control; A Func- 
t ion Space Approach , 19. Recknan, G. T,oos , and 
H.p. kun/i , eds., Spr 1 nger-Verl ag , Berlin, 1973. 
-lecture Notes in Economics and Mathematical 
Systems, vo! . 84.1 

: ^Ba lakrishnan, A.V., “Modelling and Iden- 
tification Theory: A Flight Control Applica- 

r - ! 0 n * ” _ a n(1 Appl icdti ons of Variable 

Struc ture" Sy stems, R.B. Mo;>ler and A. Robert i , 
e<1s.. Academic Press, New York, 1972. 

^^Iltff, Kenneth W., " Ident i f i cat ion and 
Stochastic Control with Application to Flight 
Control in Turbulence, 1 ' Ph.D. Di ssertat i on , Univ. 
of California, Los Angeles, May 1973. 

^^Maine, Richard E., and Iliff, Kenneth W., 

” The Theory and Practice of Estimating the 
Accuracy of Dynamic FI i ght-Determi ned Coef- 
ficients," NASA RP-1077, 1981. 

^"Parameter Estimation Techniques and 
Applications in Aircraft Flight Testing," NASA 
TN D- 764 7, 1974. 

^"Methods for Aircraft State and Parameter 
Identification," AGARD-CP-1 72, May 1976. 

^ 'Weuhaus , W.U., "Summary of Dynamic 
Stability and Control Flight Research Conducted 
Utilizing a B-25J Airplane," Rep. no. TB-405-F-K1, 
Cornell Aeronautical Laboratory, Buffalo, New 
York, May 1948. 

^Seamans, R.C., Jr., Blasingame, B.P., and 
Clement son, G.C., "The Pulse Method for the 
Determination of Aircraft Dynamic Performance ," 

J. Aeronaut. Sc i . , vol . 17, no. 1, Jan. I960, 
pp. 22-38. 

^Greenberg , Harry, "A Survey of Methods for 
Determining Stability Parameters of an Airplane 
From Dynamic Flight Measurements," NACA TN-2340 
1961. 

^Shinbrot, Marvin, "On the Analysis of Linear 
and Nonlinear Dynamical Systems From Transient- 
Response Data," NACA TN-3288, 1954. 

^Howard, J., "The Determination of Lateral 
Stability and Control Derivatives From Flight 
Data," Can, Aeronautics Space J., vol. 13, no. 3, 
Mar. 1967, pp . 126-134. 

^Wolowicz, Chester H., "Considerations in the 
Determi nat i on of Stability and Cont r ol Derivatives 
and Dynamic Cha rac ter i st ics From Flight Data," 
AGARD-AR-649, Part 1, 1966. 


33 Rampy, John M., and Berry, Donald T., 

Determi nat i on of Stability Derivatives From 
Flight Test Data by Means of High Speed Repetitive 
Operation Analog Matching," FTC -TDR -64-8, Edwards 
Cal i fornia , May 1964. 

34 Rynaski , Edmund G. , "Application of Advanced 
Identification Techniques to Nonlinear Equations 
of Motion, Proc. Stal 1 /Post-Sta 1 1 / Sp i n Symposium, 
Wright-Patterson Air Force Base, Ohio, Dec. 15-17 
1971, pp. 0-1 to 0-18. 

3 ^Iliff, Kenneth W., and Taylor, Lawrence W., 
Jr., "Determination of Stability Derivatives from 
Flight Data Using a Newton-Raphson Minimization 
Technique," NASA TN D-6579, 1972. 

36 Taylor, Lawrence W., Jr., and Iliff, Kenneth 
W., "A Modified Newton-Raphson Method for Deter- 
mining Stability Derivatives From Flight Data," 
Second International Conference on Computing 
Methods in Optimization Problems, San Remo, Italy, 
Sept. 9-13, 1968, Academic Press, New York, 1969, 
pp. 353-364. 

37 Larson, Duane B., and Fleck, John T., 
"Identification of Parameters by the Method of 
Quasi 1 inearization," Report 164, Cornell Aero- 
nautical Laboratory, Buffalo, New York, May 1968. 

3 ^Bal akrishnan, A.V., Cormuni cation Theory, 
McGraw-Hill, New York, 1968. 

39 Bellman, Richard E., and Kalaba, Robert E., 
Quasi 1 inearization and Nonlinear Boundary- Va 1 ue 
Problems , American tlsevter Publishing To, Wow 

York, 1965. 

4 °Taylor, Lawrence W. f Jr., Iliff, Kenneth W., 
and Powers, Bruce, 6., "A Comparison of Newton- ~ 
Raphson and Other Methods for Determining Sta- 
bility Derivatives From Flight Data," AIAA-69-315 
Mar. 1969. ” 

41 Grove, Randall 0., Bowles, Roland L., and 
Mayhew, Stanley C., "A Procedure for Estimating 
Stability and Control Parameters From Flight Test 
Data by Using Maximum Likelihood Methods Employing 
a Real-Time Digital System," NASA TN D-6735, 1972. 

42 Ross, A. Jean, and Foster, G.W., "FORTRAN 
Programs for the Determination of Aerodynamic 
Derivatives From Transient Longitudinal or Lateral 
Responses of Aircraft," TR-75090, Royal Aircraft 
Establishment, Sept. 1975. 

43 Maine, Richard E., and Iliff, Kenneth W., "A 
FORTRAN Program for Determining Aircraft Stability 
and Control Derivatives From Flight Data," NASA TN 
D-7831 , 1975. 

44 Nagy, Christopher J., "A New Method for Test 
and Analysis of Dynamic Stability and Control," 
AFFTC-TD-75-4, Air Force Flight Test Center, 
Edwards, California, May 1976. 


fit*#' ** 


Ch f "' ff * Kenneth Maine, Richard E., and 

n»ri fe L’ Mary , F ” " Subson1c stability and Control 
Derivatives for an Unpowered, Remotely Piloted 

i&“‘ ** ™ wj^'iSir 1 ° bu '" M Fr ” ni » K 

Dot 46b1 '''^ ~ ,e * G ” " A Corre1a t>on Between Flight- 

he P r?Sr B Va !^ es and w ’ od-Tunnel Data for 
the X-24M Research Aircraft," NASA SX-3371 , 1976. 

47 Deylum Paul M., "AFFTC Experience in the 
Use >t Flight-Test Derived Stability Derivatives " 

• aper 14, AGARD Fluid Dynamics Panel Symposium on 

iayT-?;; 3 !^ Paraneters * Athe "*’ Greece, 

4H Hollenan, Euclid C., "Summary of Flight 
lull De , tern1n ? Spin and Controllability 
n/H? t K. St ’ C ! ° f 3 * en0tel * P ’ loted > Large-Scale 
1976 F,yhter A, rplane Model," NASA TN D-8052, 

4y Smith. Harriet J., "Flight-Determined 
Stability and Control Derivatives for an Execu- 
tive Jet Transport," NASA TM X-56034, 1975. 

.. ^° Su1 t, William T., "Aerodynamic Parameters of 

NASA^TN* 0-6643^' 1972^ X ^ raCteC * Fr0 " F1 Data -“ 

Slfrp^ D.R., "Practical Applications of 
Parameter Identification," A1AA-77-11 36, Aug. 

v Ross, A. Jean, "Determination of Aerodynamic 

muht at " V Methor I ran *! er,t Response in Manoeuvring 
' 9 ?* ^Lhods for Aircraft State and Parameter 

to 14 10 AGARD-CP-1 72, May 1975, pp. 14-1 

“Gould, o,g., and Hindson, W.S., "Estimates 

° v/cim ''xty Derivatives of a Helicopter and 

Ai re raft Si'?? FHght [,ata *" Methods for 

btate and Parameter Identification 

GARO-CP-1 72, May 1975, pp. 23-1 to 23-9. 

Kliln, V., Longitudinal Aerodynamic Deriva- 

EltHrlLV 16 ?? 0elta ' Win 9 Research Aircraft 
Extracted From Flight Data," CIT-FI-74-023 

Cranfield Inst, of Technology, July 1974 . 

SSschuetz, A.J., “Low Angl e-of-Attack 

Aerodynamic Parameters of Navy T-2 
framer Aircraft Extracted From Flight Data: A 

Comparison of Identification Technique. Volume I 
7 Data Acquisition and Modified Newton-Raphson 
Analysis " NADC-74181-30-V0L-1 AD-A013181. Naval 
«1une^976? Pment Center ’ Har *’"ster, Pennsylvania, 

56 Marchand, M., and Koehler, R., Determina- 
tion of Aircraft Derivatives by Automatic Param- 

Method^ J f Stl, T t an r Freq,ienc >' Response Methods," 
•ethoris for Aircraft State and Parameter Ident- 
•fication, AGARD-CP-1 72, May 1975, op. 17-f t0 

Dynamic Stability Parameters," A6ARD- 
CP-235, 1978. 


1979 58 Parameter Ident1ficat1on ." AGARD-LS-104, 

tifi!!tS?of R r an V “ Maximu,n Likelihood Iden- 
ill? 5 * of Aircraft Parameters," Eleventh 
.o nt Automatic Control Conference of the Ameriran 

s ? ntro1 r »'' im. 

M . 6 ° T ?®r, James S., Powell, J. David, and 

Advanced^erh*’ ’ " Th f Use of Snooth ’ n 9 and Other 
Advanced Techniques for VTOL Aircraft Parameter 

Identification. Final Report,” Naval Air Sucto»c 
Command Contract No. N00019-69-C-0534 Svst^I 
Control , In.., Cali fornia, June' 1970. 

hiiit lG ^r)ach t O.H., "The Determination of Sta- 

itrtirc eriv ^ t ^ ve ? and Performance Character- 
tics From Dynamic Manoeuvres R pd VTH iai 

Delft Univ. of Technology, Dept. of Aerosolce’ 
Engineering, Delft, The^theMands, ^ 1971 . 

honkers, H.L., "Application of the Kalman 

Tes^Datl Reconstruction From Flight 

SaS Fvlnr 1 ? 1 9 Est1mation of Instrumental 
Bias Error Corrections," Report VTH-162 Delft 

univ. of Technology , Dept, of AerospaceEn*!- 
neering. Delft, The Netherlands, Feb. 1976. 

63 Chen, Robert T.N., and Eulrich. Bernard J 
Parameter and Model Identification of nu-i • J ” 
Oyna.l.,1 Systems UstiJ" iSSUS.1 
Smoothing Algorithm," Twelfth Joint Automatic 
Control Gonfer f nc « of the American Automatic 
PP 731-740 nC ’ 7 ' E2 ’ A “ 9 * 1971 ’ 

r enj1 -’ " 1(1e nti fication of Aircraft 
ability and Control Derivatives in tho Draco 
Of Turbulence," AIAA-77-1 1 al! Lg. w5? 

t.p Kenneth W., and Maine, Richard E., 

Estimator " °f US j n9 a Hax<n,um Likelihood 

tstimator. Methods for Aircraft State and 

£!;*!“!: Piper “• 

"F.,rth ll'lt.’ Kenneth W " and Maine, Richard E., 

E?im?? a) f e ? a ?? nS on Max1mun Likelihood 
Estimates of Stability and Control Character- 

Aug! C 1977! alned Fr ° m F1 ' 9ht Data ’” AIAA - 77 -H33, 

67 Maine, Richard E., and Iliff Kenneth u 
User's Manual for MMLE3, A General’ FORTRAN ” 
Program for Maximum Li keli hood Parameter 
Estimation," NASA TP-1563, 1980. 

“Fo™^? 6 ’ Ri ?? rd , E - a " d niff. Kenneth M.. 
Formulation and Implementation of a Practical 

Ind M 1thm for Parameter Estimation With Process 
and Measurement Noise," SIAM J. Appl . Math 
vol . 41, no. 3, Dec. 1981, pp. 558-579. ” 

b9 Maine, Richard E., and Murray, James E. 

Hn«?lM at I° n of . Pa |; an, eter Estimation to Highly 
Unstable Aircraft," AIAA-86-2020-CP, 1986/ Y 
(Also published as NASA TM-88266, 1986.) 


4 







70 Nguyen, Luat T., Gilbert William t> 

asir- 

»%^aa4s.;j:'sjriK n *^ 

E k n ”n r i:J°T eph *, W1,son - R *J*. and Enevoldson 
- . ? f Ge r r f • J o«P*. "Simulation as an Analysis 

A^£[rr pfl ->" 

River. Mar? S 16-??. I982! t,0n SW ° Slu Patuxent 

G N 73 "L n r.T n h ’ P ’ B ;’ Chi1d * «* R *. and Panageas, 
t he " Hi gMy°Mareuverabl'e^ Aircraft ° f 

snsaris?; Research 

gtli^rS’HSly ia‘ne^ra^e^>c7aft SyStemS 

Technolog: (Hi NAT) Vehicle." AIAA-79-1 789 ,^ 


II!! 1 Pert ’ nent values as a function of 

L Lp(L) 


0 

1 

2 

3 

4 


-0.5000 

-0.3005 

-0.2475 

-0.2500 

-0.2500 


L4(L) 

15.00 
9.888 
9.996 

10.00 
10.00 


♦a> 

0.9048 

0.9417 

0.9517 

0.9512 

0.9512 


*(L) 


Jl 


2.855 

1.919 

1.951 

1.951 

1.951 


21.21 

0.5191 

5.083 x 10-4 
1.540 x 10-9 
1.060 x 10-14 


“HiMAT^erodynamk ^esign^n^n ?ght r9e 

Experience," AIAA- 8 ! -2 4 !3, llSJ . ^ 

and Control ’Flight a Test’Terhn Ce Shu !i t,e Stability 
Aug. 1980. 9 1 Technit loes,“ AIAA-80-1608. 

"The fc“"wro% C terM!n»?rr M- J '” y 

»««,? •sss'i sssis# 

ii. ®We i 1 , Joseph, and Powers, Bruce G 
Part 

Configurations/^^SA^TM-SlBei ! & 198K *" 9 

Table 2 Pertinent values as 
a function of iteration 


0 

1 

2 

3 

4 
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Lp(L } 

L«(l) 

♦(L) 

*(U 

Jl 

-0.5000 

-0.3842 

-0.3518 

-0.3543 

-0.3542 

-0.3542 

15.00 

10.16 

10.23 
10.25 

10.24 
10.24 

0.9048 

0.9260 

0.9321 

0.9316 

0.9316 

0.9316 

2.855 

1.956 

1.976 

1.978 

1.978 

1.978 

30.22 

3.497 

3.316 

3.316 

3.316 

3.316 


Jable 3 Mean and standard deviations for estimates 
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ulL P ) deviation, 0( L p ; 


5 

10 
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Standard deviation 
of the sampl e mean, 
o(Lp)/VfI 
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0.0336 

0.0196 

0.0129 
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Cramfr-Rao bound as functions of 
the square root of noise power 


Square root of Estimate 
noise power, G 0 f L p 


Cram^r-Rao 

bound 


0.0 
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0.10 

0.2 

0.4 

0.8 
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2.0 

5.0 

10.0 


*0.2500 

-0.2507 

-0.2535 

-0.2570 

-0.2641 

-0.2783 

-0.3071 

-0.3218 

-0.3975 

-0.6519 

-1.195 


0.00054 

0.00271 

0.00543 

0.0109 

0.0220 

0.0457 

0.0579 

0.1248 

0.3980 

1.279 
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